Finite-bias electronic transport of molecules in water solution 
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The effects of water wetting conditions on the transport properties of molecular nano-junctions are 
investigated theoretically by using a combination of classical molecular dynamics and first principles 
electronic transport calculations. These are at the level of the non-equilibrium Green's function 
method implemented for self-interaction corrected density functional theory. We find that water 
effectively produces electrostatic gating to the molecular junction, with a gating potential determined 
by the time-averaged water dipole field. Such a field is large for the polar benzene-dithiol molecule, 
resulting in a transmission spectrum shifted by about 0.6 eV with respect to that of the dry junction. 
The situation is drastically different for carbon nanotubes (CNTs). In fact, because of their hydro- 
phobic nature the gating is almost negligible, so that the average transmission spectrum of wet 
Au/CNT/Au junctions is essentially the same as that in dry conditions. This suggests that CNTs 
can be used as molecular interconnects also in water-wet situations, for instance as tips for scanning 
tunnel microscopy in solution or in biological sensors. 

PACS numbers; 75.47. Jn,73.40.Gk,73.20.-r 



I. INTRODUCTION 

Electron transfer in biological systems is rather differ- 
ent from electron transport in an electronic device. In ad- 
dition to the intrinsic materials differences between the 
respective conducting media, namely soft molecules in 
biology against inorganic semiconductors in electronics, 
the electron transport in these two situations differs for 
the effects and the relevance of the environment. In fact, 
while one wants to keep electronic devices in dry condi- 
tions to avoid electrostatic disorder, biological electron 
transfer is characterized by a time-evolving local dielec- 
tric environment. 

The boundary between the two fields, electronics and 
biology, becomes however blurred when one looks at the 
nanoscale. On the one hand, "conventional" nanotech- 
nology has expanded in the biology domain with a grow- 
ing number of electronic applications requiring opera- 
tion in wet conditions. In addition to scanning tunnel 
microscopy in water—, which has been available for the 
last two decades, several biological sensors have been 
recently proposed. These for instance include cancer 
markers detectors^ and protocols for DNA sequencing^. 
On the other hand, when a solution is confined at the 
nanoscale, highly ordered structures can form at room 
temperaturei. This essentially means that under such 
stringent confinement a biological system almost behaves 
like a solid. 

Interestingly, numerous experiments to date on elec- 
tronic transport through molecules are carried out in 
solution&iS. Still, with a few exceptions'"^, most of the 
theoretical calculations are performed in the dry, i.e. 
without explicitly including water molecules in the simu- 
lations. Therefore the question on how water affects the 
current- voltage characteristic in a molecular junction re- 
mains. In general one may expect that solutions made 
with polar molecules, such as H2O, may affect signif- 



icantly the transport because of the generation of local 
dipole fields. However, since the time scale of a transport 
measurement is far longer than the typical molecular re- 
arrangement of a solution, one should ask what is the 
time-averaged dipolar field near the molecule of interest. 
This, together with the degree of localization of the rel- 
evant molecular orbitals, will determine whether or not 
the wetting conditions influence the junction electronic 
transport. 

In order to address those fundamental questions we 
have performed a number of combined molecular dynam- 
ics (MD) and quantum transport simulations for molecu- 
lar junctions in water. Our computational strategy is to 
investigate first the dynamics of H2O and then to eval- 
uate the electronic transport of a representative number 
of configurations, i.e. for a number of MD snapshots. In 
particular we have considered two rather different junc- 
tions, both using gold electrodes but differing for the local 
charge arrangement of the molecule of interest. The first 
is made from a polar molecule, namely benzene-dithiol 
(BDT), while the second includes locally charge neutral 
carbon nanotubes (CNTs), in particular a (3,3) metallic 
and a (8,0) insulating one. 

We find that the effect of water on the transport is 
that of effectively gating the molecule, therefore shifting 
almost rigidly its transmission spectrum. This is rather 
pronounced for BDT, but only tiny in the case of CNTs, 
and reflects the different charge distribution of the two 
classes of molecules. We associate those differences to 
the time-averaged dipole field of the water. Notably all 
the calculations are performed with density functional 
theory including appropriately corrections for the elec- 
tronic structure of water. This is an essential condition 
for quantitative predictions. 

The paper is organized as follows. In the next section 
we describe the device set-up for the calculations and 
briefly the computational tools used. Then we analyze 
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the transport. First we look at the electronic structure 
of a gold capacitor with water wetting the two electrodes 
and then we consider finite bias conductance across BDT 
and the CNTs. Finally we present our main conclusions. 



II. METHODS 

Classical MD calculations are performed with the 
NAMd2 packag©i2. BDT is attached to the two gold elec- 
trodes at the Au(lll) hollow site, which has been pre- 
viously calculated to be the low energy bonding position 
for BDT on Au(lll)JJ'. We define our coordinate system 
with the z axis along the (111) direction and the x-y plane 
orthogonal to it. All the calculations are carried out with 
periodic boundary conditions in the x-y plane. TIP wa- 
ter molecules are added to a (25 x 17.3 x 7.3) A^ box 
intercalated between the Au electrodes. Note that the 
dimension along the z axis is chosen in such a way that 
the Van der Waals distance between the Au(lll) surface 
and the water molecules is taken into account. In total 
there are 83 II2O molecules in the simulation box. Since 
our main objective is that of examining the effects of wa- 
ter over the conductance of a Au/molecule/Au device, 
we always fix the atomic positions of both the molecule 
and the electrodes. 

The interaction between the H2O molecules and gold 
is treated at the level of a 12-6 Lcnnard- Jones potential 




with parameters for Au (e = 0.039 kcal/mol and a — 
2.934 A) taken from the literature [12]. Periodic bound- 
ary conditions are applied with a cutoff of 12 A for 
long-range interactions. In order to maintain the size 
and shape of the cell constant, we perform simulations 
in the micro-canonical ensemble, with re-initialized ve- 
locities to 300 K for every 1000 time-steps and with 
a time-step of 2 fs. The trajectory is recorded ev- 
ery 4 ps from the initial equilibration of 1 ns to a to- 
tal simulation time of 20 ns. For the systems com- 
prising the CNTs, Au(lll)/CNT(3,3)-H20/Au(lll) and 
Au(lll)/CNT(8,0)-H2O/Au(lll), and for the reference 
Au(lll)/H20/Au(lll) capacitor, similar conditions and 
procedures are followed. The parameters for the aromatic 
carbon atoms are used for CNTs. 

For each system we calculate the transport proper- 
ties for a set of representative MD configurations taken 
after equilibration. From these we can then estimate 
the fluctuations in the transmission and current over 
time, as well as their time-averaged values. We use the 
SMEAGOL ah initio electronic transport code to calculate 
the zero bias transmission coefficients and the current- 
voltage iJ-V) characteristics. SMEAGOL combines the 
non-equilibrium Green's function method with density 
functional theory (DFT)^^"— and has the pseudopoten- 
tial code SIESTA^^ as its electronic structure platform. 



The local density approximation (LDA) to the ex- 
change and correlation functional is adopted throughout. 
The atomic self-interaction correction (ASIC)^^ scheme 
however is used for the water and the BDT molecule, 
in order to bring their ionization potentials (IPs) in 
closer agreement to experiments. This has been already 
proved to be a successful strategy for aligning correctly 
the highest occupied molecular orbital (HOMO) energy, 
fHOMO I of the molecule under consideration to the Fermi 
level (i^p) of the electrodesiiii^. Note that the same 
corrections also reproduce well the band-gap of many 
insulating oxides, after the ASIC potential is rescaled 
appropriately^^. Such a rescaling in bulk crystals is at- 
tributed to charge screening, which in solids is usually 
stronger than in molecules. In general we use a scal- 
ing parameter, a, ranging between and 1, to adjust 
the amount of self-interaction correction included (a = 
corresponds to the LDA, a = 1 is the full ASIC). Usually 
a is 1 for small molecules, it is around 0.5 for insulating 
oxides and it vanishes for metals. For this reason ASIC 
is never applied to Au. 

For the transport calculations we use a real space mesh 
cutoff of 200 Ry and an electronic temperature of 300 K. 
The unit cell includes 5 Au atomic layers on each side of 
the Au(lll) surface, which are enough to screen charg- 
ing at the Au-molecule interface. In order to reduce the 
system size the Au 5d shell is kept in the core, so that 
we consider only 6s orbitals in the valence. We use a 
single-f^ basis for Au 6s, specifically optimized to give 
the correct work function for the Au(lll) surface. The 
rest of the basis set is double-^ for the C s and p orbitals, 
and double-^ plus polarization for S (s and p). For the 
H2O molecules we use a double-C basis for both H and 
O. However when calculating transport through CNTs 
we reduce it to single-C, in order to keep the size of the 
density matrix tractable. 

The charge density is obtained by splitting the integral 
of the Green's function into a contribution calculated 
over the complex energy plane and one along the real 
axiai^ii^. The complex part of the integral is computed 
by using 16 energy points on the complex semi-circle, 16 
points along the line parallel to the real axis and 16 poles. 
The integral over real energies necessary at finite bias is 
evaluated over at least 1000 point o^^'^'^ per eV. 



III. RESULTS AND DISCUSSION 

A. Au capacitor in water 

In a transport calculation it is crucial to describe ac- 
curately both the electrodes' work function, W ^ and the 
IP of the molecule, so that the correct alignment of 
the molecular levels to the electrodes' E-p is reproduced. 
With this in mind we first analyze the electronic struc- 
ture of water sandwiched between Au electrodes. Such a 
setup corresponds to calculating the electronic structure 
of a Au parallel plate capacitor, where the two plates 
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FIG. 1: (Color online) An capacitor with water as dielectric 
medium. The top panel shows the unit cell used for the MD 
simulations, which includes H2O confined between the two 
Au electrodes. In panel (a) we show the probability to find 
a O (solid curve) or a H (dashed curve) atom at a given z 
position. Panels (b-d) are the planar averages of the Hartree 
electrostatic potential as function of position along the 
transport direction, z, for different setups: (b) entire junction, 
(c) junction without the electrodes, and (d) junction without 
water. The grey curves, merging in a shadow, are the results 
for 21 snapshots taken at different times after equilibration, 
the black solid curves are their time averages. The dashed 
curve in (c) shows the difference between the time-averages 
of panels (b) and (d) . 



are separated by water. The unit cell, containing 408 
H2O molecules and 480 Au atoms, is shown in Fig. [1] 
The statistical distribution of H2O about the Au plates 
is described in Fig. [ija), where we plot the normalized 
probability, ^(z), to find O (solid curve) or H (dashed 
curve) at a given position z in the cell (/ p(z) dz = 1). 
Such a distribution is obtained by using the MD data for 
all the time-steps included in the 16 ns simulation after 
equilibration. We note that p(z) is constant, in the mid- 
dle of the gap between the plates, indicating an average 
random arrangement of the H2O molecules. In contrast 
close to the Au interface there are marked oscillations in 
p(z), signaling a correlation of the water position with 
respect to the Au surface. Note that the peaks in p(z) 
are found at the same positions for O and H atoms. This 
indicates that on average there is no net dipole at the 
interface. 

Next the Au work function is calculated by using the 
same Au parallel plate capacitor of Fig. [T] after we have 
removed the water, i.e. with vacuum as spacer between 
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FIG. 2: (Color online) DOS projected onto II2O for the capac- 
itor of Fig. [1] and obtained for different values of the ASIC 
scaling parameter a. The grey lines are the superimposed 
curves for 21 MD snapshots and the solid black one is their 
average. The dashed curve is the average DOS for the same 
MD snapshots for a system where Au is replaced by vacuum. 

the plates. W is then the energy difference between 
and the vacuum potential, Vvacuum- Such an exercise is 
reported in Fig. [Ijd) where we show the planar average 
Vh of the electrostatic Hartree potential along the x-y 
plane. Note that the absolute value of Vn in Au is arbi- 
trary, since it depends on the pseudopotentials. However, 
in the middle of the capacitor Vh = Vvacuum- There- 
fore, setting Kacuum = we have W — —Ep. We find 
W = 5.3 eV, in good agreement with experiments. 

We now turn our attention to the electronic structure 
of water. This must be extracted for its liquid phase, 
i.e. from the MD simulations. We consider data aver- 
aged over 21 structural configurations corresponding to 
21 equally spaced MD snapshots. The so-calculated pla- 
nar average of the electrostatic potential is shown in Fig. 
Wljo). The grey lines, merging in a shadow, are the super- 
imposed curves for all the snapshots, while the black line 
is their time-average. We note oscillations of t^j close 
to Au. These are due to the arrangement of the water 
molecules with respect to the Au surface. In contrast 
in the middle of the junction Vh is rather flat due to 
the average random orientation of the molecules [see also 
Fig. Ha)]. 

From the same simulations we can extract the IP of liq- 
uid water. In transport one wants chomo of the molecule 
of interest to correspond to the actual negative of its 
jp20,2i^ Although this should be the case for exact DFT, 
it happens only rather rarely in practice for the stan- 
dard local approximations of the exchange and correla- 
tion functionaP^. In the case of liquid H2O the experi- 
mental value of the IP is in the range between 9.3 and 
10 eV, in rather good agreement with recent ab initio 
calculations of the electron removal energy (see [23,23] 
and references therein). 

We evaluate chomo for liquid H2O by calculating 
the density of states (DOS) projected onto the water 
molecules. This is shown in Fig. [5] (Kacuum = 0), where 
again the grey curves are the superimposed data for all 
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the MD snapshots and the black one is their average. We 
note that there are large fluctuations in the DOS over 
time with energy shifts of the order of 1 eV. The aver- 
age value however is smooth and can be reliably taken as 
the liquid water DOS. In this way the top of the water 
"valence band" is calculated to be only -6 eV in LDA, 
in agreement with previous DFT calculations^^, but still 
far from the experimental value. 

We then apply ASIC to H2O and find that ehomo 
moves to -8.5 eV for a = 0.5 and to -11 eV for a = 1.0, 
so that a — 0.7 fits the average experimental value (- 
9.5 eV). Such an optimal value, as usually with ASIC—, 
in general improves the entire electronic structure and 
returns an HOMO-LUMO gap of about 6.4 eV, in good 
agreement with the experimental value of 6.9 eV.-^S Note 
that a = 0.7 is typical for moderately ionic insulating 
oxidesii. 

In order to analyze the effects of the AU/II2O inter- 
face over the water IP and DOS we perform a second set 
of calculations, where we remove the Au plates and we 
replace them by vacuum. This corresponds to an hypo- 
thetical H2O slab. In this case [see Fig. [T{c)] Vh at a 
given MD time-step has a finite slope in the vacuum re- 
gion, which is caused by the non-compensated dipoles at 
the H2O external surface. These dipoles produce a long- 
range electric field outside slab, so that Vvacuum of a sin- 
gle snapshot is not defined. However the time-averaged 
Vn (black curve) is approximately flat in the vacuum, 
demonstrating that, although at each time-step surface 
charge may lead to long-ranged electric fields, its time- 
average is actually zero. 

If we now take the average Vn away from the water 
molecules as Vvacuum, we can plot the time-averaged DOS 
for the water slab and superimpose it to that calculated 
for the AU/H2O/AU capacitor [see Fig. [2]. We find that 
the two DOSs overlap on each other, confirming our re- 
sults for the water IP and the fact that on average there 
is little electronic interaction between Au and H2O. Our 
results also suggest that one should ideally use periodic 
boundary conditions to simulate the electronic structure 
of molecules in solution. These eliminate the possible 
spurious electric fields in the vacuum, which can lead to 
an unphysical rearrangement of the energy levels. Fur- 
thermore for simulations of H2O surfaces it is essential 
to consider time-averages, so that the water electric field 
vanishes in vacuum. Figlljc) also shows the difference 
between the time-averaged Vn of the Au capacitors with 
and without water (dashed line). This difference is al- 
most identical to the time-averaged Vh for the water slab 
and is consistent with the fact that the Au electrodes do 
not induce any noticeable change in the average DOS of 
H2O. 



IV. BENZENE-DITHIOL 

Electron transport through BDT attached to 
Au electrodes has been extensively studied both 




FIG. 3: (Color online) Unit cell used for the BDT molecule 
in water attached to Au electrodes. 

experimentally^^"— and theoretically . ^ ^1 In fact, 
because of its simple structure, BDT is an ideal system 
for comparing theory with experiments. However, also 
for BDT the LDA description is not adequate and it is 
necessary to employ ASIGii^i^. In order to limit the 
number of adjustable parameters we set the same a 
for both BDT and H2O, and check that such a value 
reproduces well previous transport calculations for the 
Au/BDT/Au junction in dry conditionsiiii^. 

We first investigate the V = transport (the cell used 
is shown in Fig. [S]). In Fig. |4]the transmission coefficient 
T{E; V — 0) for one MD snapshot is shown as function 
of energy, E, for different values of the ASIC scaling pa- 
rameter a (solid curves) . The same quantity is compared 
to that calculated for the same cell, this time without in- 
cluding water (dashed curves). In general the effects of 
water are two-fold: firstly there is a shift of the BDT 
transmission peaks to lower energies, and secondly there 
appear additional sharp transmission peaks, which are 
attributed to resonant transport through the electronic 
states of H2O. Interestingly the height, the width and the 
relative position of the transport peaks with respect to 
each other is unchanged when water is present. 

Therefore the main effect of adding water is to shift 
the energy levels of BDT, so that water acts as an exter- 
nal gate. Since the BDT molecular orbitals extend over 
the entire molecule and are strongly coupled to Au^^, all 
the levels shift by approximately the same amount. If 
the levels were more localized, we might have expected 
a change in their relative position, sensitively dependent 
on the local configuration of H20'^°. 

We now analyze in more detail the electronic proper- 
ties of H2O for the same MD snapshot. In Fig. [S] the 
DOS projected onto the water molecules is shown over 
the same energy range as that of T in Fig. |4l It is clear 
that the additional peaks in the transmission (Fig. S]) are 
at energies where H2O has a finite DOS, confirming that 
these are due to transport through the electronic states 
of water. The additional transmission peaks below E-p 
are rather close to the Fermi level when a = 0, whereas 
they move down in energy as a is increased. In contrast 
the position of the peaks above Ep is almost constant 
for different a, since the ASIC mainly affects occupied 
states. We also find that for a = the water HOMO 
is about -6.3 eV from Kacuum, while it is at -8.7 eV for 
a = 0.5. These values agree with the values obtained for 
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FIG. 4: (Color online) Transmission coefficient for transport 
through BDT in water as function of energy for one MD snap- 
shot and for different values of the ASIC scaling parameter a. 
The dashed curves are for BDT in dry conditions (no water), 
while the solid curves are for BDT in solution. 
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FIG. 5: Density of states projected onto II2O for the junction 
of Fig. [3] as function of energy for one MD snapshot and for 
different values of the ASIC scaling parameter a. 
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FIG. 6: (Color online) Current, I, as function of the bias volt- 
age, V, through the BDT molecule calculated at the first MD 
time-step and for different values of the SIC scaling parameter 
a. The dashed curves are for the molecule in the dry while 
the solid curves are for BDT in solution. 
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FIG. 7: (Color online) Time averaged (a) transmission coeffi- 
cient and (b) DOS projected onto the H2O molecules for BDT 
attached to Au. The calculations are obtained with ASIC and 
a = 0.7; the solid black curves are the time-averages over 201 
time-steps while the grey ones forming a shadow are for each 
of the 201 time-steps. The dashed curve is for the dry situa- 
tion. 



the water slab. 

The self-consistent current-voltage, I-V, curve is 
shown next in Fig. [5] for the same MD snapshot and for 
different values of a. Generally speaking the presence of 
water leads to a reduction of the current, which is more 
pronounced for small a. The reason for such a reduction 
is that ehomo is very close to E-p for small a, so that a 
tiny shift of the BDT levels to lower energies consider- 
ably reduces T{Ep;V « 0). For a = 1 there is almost 
no change in the current, since E-p is approximately at 
mid-gap already in dry conditions. 

We now briefly discuss the most appropriate value of 
a in this context. For BDT in the gas phase a = 1 gives 
ehomo close to the experimental -TPj^^d^ However, when 
BDT is immersed in water additional screening lowers 
down the value of a. We then take a — 0.7 (the optimal 
value for H2O), which provides a good IP for water and 
also accounts for the additional screening in BDT due to 
the solution. Importantly our results depend little on the 
exact choice of a, as long as it is of the order of 0.5, i.e. 
such that ehomo for H2O is well below the Au Ep. Note 



that using a = would lead to the erroneous prediction 
that water becomes conducting at about 1 V of bias. 

We now move to calculate the time-averaged trans- 
mission coefficient and the I-V curves. In this case 
T{E; V = 0) and the H2O DOS are evaluated over 201 
snapshots taken in the last 16 ns of our MD simulations, 
while the I-V curves are evaluated over only 21. Trans- 
mission and DOS are presented in Fig. [71 where again 
the curves for the single snapshot calculations are plot- 
ted in grey to form a shadow, while their average is a 
solid black line. In general, when H2O is introduced in 
the simulation, there is a rigid shift of the entire spec- 
trum towards lower energies with respect to the dry sit- 
uation (dashed line). This is because BDT and H2O are 
both polar molecules and in time the water molecules 
arrange around BDT so to screen the local dipole field. 
Such screening moves the average BDT molecular levels 
to lower energies. 

This analysis is confirmed in Fig. [SJ where we present 
the O and H position distributions along the y direction, 
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FIG. 8: (Color online) Probability to find an O (solid curve) 
or a H (dashed curve) atom at a given position y in the 
Au/BDT/Au junction, for atoms whose x coordinate lies 
within the shadowed region. In the top panel we show a rep- 
resentative snapshot of the atomic configuration of the water. 
Note that in the first solvation layer the H2O molecules are 
oriented with the H atoms pointing towards the BDT. 



for those H and O atoms lying either above or below the 
plane of the BDT (shadowed region in Fig. [8|). Note that 
we define the y axis as the direction perpendicular to the 
plane of the BDT. In contrast to the case of the Au ca- 
pacitor, now the p(2/)'s of O and H ions differ near BDT. 
In particular we find that H approaches BDT approxi- 
mately 1 A closer than O. This means that on average 
the first solvation layer is oriented with the H atoms of 
H2O molecules pointing towards the BDT, as suggested 
by elementary electrostatics, since the C and S atoms are 
fractionally negatively charged, while the H atoms have 
a positive charge. The second peak of the H atoms over- 
laps with the first peak of O, indicating that while one 
of the H atoms points towards the BDT, the second one 
aligns with the negative O atoms. It is also interesting 
to note that the O distribution has a second pronounced 
peak in addition to that close to the BDT, signaling a 
relative large degree of order also in the second solvation 
layer. 

We finally turn our attention to the fluctuations. Gen- 
erally, time-fluctuations in the position of the BDT 
single particle levels result in zero-bias conductance 
fluctuations^. These cause both a reduction in the av- 
erage height of the various transmission peaks and at the 
same time an increase of their widths. For BDT attached 
to Au this second effect is rather small, since the trans- 
mission peaks at each time-step are already rather broad, 
due to the strong electronic coupling with the electrodes. 
However we expect that for molecules weakly coupled to 
the electrodes and thus presenting sharp peaks in T{E) 
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FIG. 9: (Color online) Histogram of enoMO extracted from 
the maximum in T{E) at around -2 eV (see Fig.[7| for a — 0.7 
and for 201 MD time-steps. A'' is the number of times enoMO 
is found in the given energy window specified by the bin width. 
The solid line indicates the time-averaged enoMO, while the 
dashed red one marks the result in dry conditions. 
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FIG. 10: (Color online) I-V curve for a BDT molecule at- 
tached to gold in presence of water obtained with ASIC 
{a — 0.7); the solid black curve corresponds to the aver- 
age current over 21 time-steps, the grey curves merging in 
a shadow are the I-V curves of each individual configuration 
and the red dashed curve is for BDT in dry conditions. 



this effect will be more pronounced, probably dominating 
the energy level broadening. 

As already mentioned before the electrostatic interac- 
tions of water with the BDT mimics a gate potential. At 
each MD time-step such an effective gate voltage changes, 
depending on the relative position of H2O. In order to 
quantify the fluctuations of the BDT molecular levels we 
track the position of ehomo [from the peak in T{E)] over 
time, and display the result in the form of a histogram 
in Fig. [9] In the plot N is the number of counts ehomo 
is found in a particular energy window, the dashed red 
line indicates the position of ehomo in the dry, while 
the solid black line indicates the time averaged enoMO in 
H2O solutions. Clearly ehomo fluctuates between -1.8 eV 
and -2.4 eV, i.e. in an energy range of 0.6 eV. The time- 
averaged enoMO is about 0.6 eV below ehomo for the dry 
molecule, which means that the effective water-induced 
gating potential is about 0.6 eV. 

Finally in Fig. [TU] we present the self-consistent I-V 
characteristics, where we conclude that the current fluc- 
tuates about its time average by approximately ±20%, 
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while it is reduced from that in the dry by about 35%. 
This discussion is based on ASIC calculations with the 
optimal value of a = 0.7. As discussed previously, the 
most appropriate correction for BDT in the dry is a ~ 1. 
If the same correction is exported into the wet situation 
the reduction of the current due to water becomes almost 
negligible. 



V. CARBON NANOTUBES 

We now move to the analysis of the effects of II2O- 
wetting on non-polar molecules, i.e. molecules present- 
ing local charge neutrality. In particular we choose two 
different CNTs: i) (3,3) metallic armchair and ii) (8,0) 
insulating zigzag. 

In the case of the metallic (3,3) CNT we again per- 
form 20 ns long MD simulations and calculate the observ- 
ables over 201 equally spaced snapshots in the last 16 ns. 
The unit cell used is shown in Fig. 1111 for one particular 
MD snapshot. This has a (20.0x17.3) cross section 
and contains 480 Au atoms, 192 C atoms and 360 II2O 
molecules. The Au-CNT distance is simply obtained by 
adding the Au and C atomic radii (respectively 1.44 A 
and 0.7 A) and it is close to that obtained by total energy 
minimization^. We note that the exact conformation of 
the Au-CNT bonding is not known and that changes in 
bond structure lead to quantitative changes in the trans- 
mission spectral^ Here however we are mainly interested 
in investigating how the transmission is affected by the 
water so that the precise bonding geometry is less impor- 
tant. 

As already mentioned before, because of the large sys- 
tem size here we use a single-C basis for II2O. We 
verified that this gives a similar IP to that obtained 
with the double-C basis. In what follows we will use 
a = 0.7 for H2O, but no ASIC for the CNTs, since 
their electron screening is good. We have verified that 
the band-structure of the (3,3) CNT agrees well with 
previous calculations 1^ In particular we obtain a CNT 
work function of 4.4 eV in good agreement with previous 
calculations2i. The IP for (3,3) CNT is not available ex- 
perimentally, but that of similar CNTs ranges between 
4.8 eV and 5.0 eV^^— thus is not far from what calcu- 
lated here. 

Since W of Au is about 1 eV larger than that of the 
CNT, electrons transfer from the CNT into Au, leading 
to a substantial band-bending. This is demonstrated in 
Fig.lTlTa). where the planar average of the Hartree poten- 
tial is plotted for the Au/CNT/Au junction in dry condi- 
tions. Close to the Au/CNT interface the oscillating Vh 
is higher than in the middle of the junction, whereas for 
an infinite CNT Vh oscillates around a constant average. 
We note that such charging effects have been neglected in 
previous tight-binding calculations. 3^ Importantly how- 
ever charging leads to a shift in the transmission spec- 
trum, so that it is important to include such an effect in 
a self-consistent way. 
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FIG. 11: (Color online) Au/CNT/Au junction. The top panel 
shows the unit cell used for the MD simulations, which in- 
cludes the CNT, Au electrodes and water molecules. In pan- 
els (a) and (b) we present the planar averages of the Hartree 
electrostatic potential, Vh, as function of position along the 
transport direction z: (a) junction without water, and (b) 
junction with water. The grey curves, merging in a shadow, 
are the results for all 201 MD snapshots, the black solid curves 
are their time-averages. 



Next we look at the wet situation of Fig. ITlT b). In 
this case Vh for a single MD snapshot oscillates dramati- 
cally along the CNT. However, when the time average is 
considered [solid black curve in Fig. [TlTb)] a regular pat- 
tern emerges, where Vh resembles closely that obtained 
in the dry. This confirms that on average the position 
of the H2O molecules away from the interface is random. 
Since CNTs are hydro-phobic, we expect the interaction 
between the water and the CNT to be weak. This is con- 
firmed by taking the difference between Vh calculated 
with and without H2O molecule and observing that the 
resulting curve matches closely that of the water slab 
calculated previously [see Fig. [Ijc)]. 

The transmission coefficients for all the 201 MD snap- 
shots are shown in Fig. [12] as super-imposed grey curves 
together with their average (solid black line), the same 
quantity calculated in dry conditions (dashed red line) 
and the total number of open channel in the CNT (dot- 
dashed green line). In the figure we shift E-p in such 
a way that Vh for the infinite CNT matches Vh of the 
CNT attached to Au without water in the middle of the 
junction. The necessary shift is about 0.8 eV, which cor- 
rectly corresponds to the difference in the work functions 
between the CNT and Au. 

The main result is that the average transmission in 
wet conditions and that of the dry junction overlap al- 
most exactly, demonstrating that in this case H2O has 
no gating effect. This can be easily understood by recall- 
ing that, since the CNT has no polar edges, the average 
H2O conformation presents no net electrical dipole, so 
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FIG. 12: (Color online) Transmission coefficient as function 
of energy for the Au(lll)/CNT(3,3)-H20/Au(lll) junction. 
The solid black curve corresponds to the average T over 201 
MD snapshots, the dashed red curve is the transmission for 
the same junction in the dry and the dot-dashed green line 
indicates the number of channels (spin-degenerate) for the 
infinite CNT. The grey curves, merging in a shadow, are 
T{E; V = 0) for each of the MD snapshot. 

that on average there is no shift of the CNT energy lev- 
els. Of course, each individual MD snapshot displays a 
dipole and the CNT energy levels get shifted. This leads 
to fluctuations in the transmission. As a result of the 
dipole fluctuations, we find that that some of the sharp 
transmission peaks visible in the dry are broadened up 
and have an average reduced height in solution. In some 
extreme cases (see for instance the sharp peak at about 
-1 eV) they are completely washed out by the fluctua- 
tions. 

Again in order to quantify the fluctuations of T{E), we 
choose a particular molecular level (transmission peak) 
and follow its time fluctuations. Here we select the well- 
defined peak at -1.7 eV below Ep and present its energy 
distribution histogram in Fig. 1131 This time the peak 
position fiuctuates over the tiny energy range of 0.06 eV, 
which is one order of magnitude smaller than that of 
the HOMO of BDT (see Fig. E]). The origin of such 
small fluctuations is twofold: firstly the interaction be- 
tween H2O and the CNT is very weak due to the hydro- 
phobic nature of the nanotube and secondly the CNT 
TT-like molecular states are delocalized, so that local fluc- 
tuations in the electrostatic potential largely cancel out 
over the entire molecule. We also find that the differ- 
ence between the average peak position (solid black line 
in Fig. [13]) and that in the dry (dashed red line) is only 
0.01 eV. This is also much smaller then the same quantity 
calculated for BDT (0.6 eV). 

Finally we discuss results for the insulating (8,0) CNT. 
The simulation cell is identical to that of the (3,3) 
case, but this time we have 288 C atoms and 322 H2O 
molecules. The MD simulations are run for 20 ns and 
only 41 snapshots are taken within the last 16 ns. We 
have reduced the number of snapshots from 201 to 41 be- 
cause this time we do not perform a detailed statistical 
analysis of the peak position. Our main results are shown 
in Fig.[T31 where we present T{E; V = 0) for all the snap- 
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FIG. 13: (Color online) Histogram of the energy position 
of the transmission peak located at -1.7 eV below Ef (see 
Fig. I12|) . The histogram has been constructed from 201 MD 
snapshots. A*' is the number of times the peak is found in the 
given energy window specified by the bin width. The solid 
line indicates the time-averaged position, while the dashed 
red one marks the result in dry conditions. 
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FIG. 14: (Color online) Transmission coefficient as function 
of energy for the Au(lll)/CNT(8,0)-H2O/Au(lll) junction. 
The solid (black) curve corresponds to the average T over 
41 MD snapshots, the dashed (red) curve is the transmission 
for the same junction in the dry and the dot-dashed (green) 
line indicates the number of channels (spin-degenerate) for 
the infinite CNT. The grey curves, merging in a shadow, are 
T{E; V = 0) for each of the MD snapshot. 



shots (grey lines), their average (solid black line), that in 
the dry conditions (dashed red line) together with the 
total number of scattering channels. The transmission 
coefficient is plotted in logarithmic scale in order to em- 
phasize the tunneling behavior in the gap. In general the 
quantitative features of the Au/CNT(8,0)-H2O/Au junc- 
tion are similar to those of the Au/CNT(3,3)-H20/Au 
one. In particular also here the average transmission al- 
most overlaps with that of the dry junction meaning that 
there is a negligible average gating. 

However for energies corresponding to the CNT gap, 
the transmission fiuctuations are rather large due to 
the tunneling nature of the transport. For instance 
T{Ef;V = 0) fiuctuates between 6.0 10"^ and 2.0 lO'^ 
within the 41 MD snapshots considered. This means that 
in tunneling conditions the presence of H2O produces 
substantial variations in the instantaneous conductance 
amplitude. However and most importantly the time- 
averaged transmission at Ep is only about 30% larger 
than that of the dry junction. This gives us the impor- 
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tant result that in general the transport in CNTs is little 
aflFected by H2O solution regardless of the metallic state 
of the CNT. 



VI. CONCLUSIONS 

In conclusion, we have investigated the effects of wa- 
ter on the transport properties of two types of molecules. 
This is done by combining classical molecular dynamics 
with ab initio electron transport calculations. Firstly, 
as an important technical result, we find that self- 
interaction corrections are fundamental for describing 
the H2O ionization energy and its band-gap. This is 
a pre-requisite for quantitative transport calculations. 
Then our main result is the finding that the H20-wctting 
conditions effectively produce electrostatic gating to the 
molecular junction, with a gating potential determined 
by the time-averaged water dipole field. Such a field is 
rather large for the polar BDT molecule, resulting in an 



average transmission spectrum shifted by about 0.6 cV 
with respect to that of the dry junction. In contrast, the 
hydro-phobic nature of the CNTs leads to almost negli- 
gible gating, so that the average transmission spectrum 
for Au/CNT/Au is essentially the same as that in dry 
conditions, regardless of the CNT metallic state. This 
suggests that CNTs can be used as molecular intercon- 
nects also in water wet situations, for instance as tips for 
scanning tunnel microscopy in solutions or in biological 
sensors. 
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